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While it is established that the effect of neutrinos on the evolution of cosmic structure is small, the upper 
limits derived from large-scale structure could help significantly to constrain the absolute scale of the neutrino 
masses. Current results from cosmology set an upper limit on the sum of the neutrino masses of ~ 1 eV, somewhat 
depending on the data sets used in the analyses and assumed priors on cosmological parameters. In this review we 
discuss the effects of neutrinos on large-scale structure which make these limits obtainable. We show the impact of 
neutrino masses on the matter power spectrum, the cosmic microwave background and the clustering amplitude. 
A summary of derived cosmological neutrino mass upper limits is given, and we discuss future methods which 
will improve the mass upper limits by an order of magnitude. 



1. INTRODUCTION 

The wealth of new data from the cosmic 
microwave background (CMB) and large-scale 
structure (LSS) in the last few years indicate that 
we live in a flat Universe where ~ 70 % of the 
mass-energy density is in the form of dark en- 
ergy, with matter making up the remaining 30 % . 
The WMAP data combined with other large-scale 
structure data |1I2| give impressive support to 
this picture. Furthermore, the baryons contribute 
only a fraction /b = flb/^m ~ 0.15 {flh and 
f2in are, respectively, the contribution of baryons 
and of all matter to the total density in units 
of the critical density Pc = 3H^/8ttG = 1.879 x 
10-29 gcm-3, where Hq = 100ft. kms'^ Mpc"^ 
is the present value of the Hubble parameter) of 
this, so that most of the matter is dark. The 
exact nature of the dark matter in the Universe 
is still unknown. Relic neutrinos are abundant in 
the Universe, and from the observations of oscilla- 
tions of solar and atmospheric neutrinos we know 
that neutrinos have a mass Si5"4'5"6'7'8''9"10' and 
will make up a fraction of the dark matter. How- 
ever, the oscillation experiments can only mea- 
sure differences in the squared masses of the neu- 
trinos, and not the absolute mass scale, so they 
cannot, at least not without extra assumptions, 
tell us how much of the dark matter is in neu- 



trinos. From general arguments on structure for- 
mation in the Universe we know that most of the 
dark matter has to be cold, i.e. non-relativistic 
when it decoupled from the thermal background. 
Neutrinos with masses on the eV scale or below 
will be a hot component of the dark matter. If 
they were the dominant dark-matter component, 
structure in the Universe would have formed first 
at large scales, and smaller structures would form 
by fragmentation (the 'top-down' scenario) . How- 
ever, the combined observational and theoreti- 
cal knowledge about large-scale structure gives 
strong evidence for the 'bottom-up' picture of 
structure formation, i.e. structure formed first 
at small scales. Hence, neutrinos cannot make up 
all of the dark matter (see e.g. JT] for a review). 
Neutrino experiments give some constraints on 
how much of the dark matter can be in the form 
of neutrinos. Studies of the energy spectrum in 
tritium decay 12 provide an upper limit on the 
effective electron neutrino mass involved in this 
process of 2.2 eV (95 % confidence limit). For 
the effective neutrino mass scale involved in neu- 
trinoless double beta decay a range 0.1-0.9 eV has 
been inferred from the claimed detection of this 
process |13ll4j . If confirmed, this result would 
not only show that neutrinos are Majorana par- 
ticles (i.e. their own antiparticles), but also that 
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the neutrino masses are in a range where they are 
potentially detectable with cosmological probes. 

The structure of this review is as follows. Sec- 
tions and 2 and 3 discuss the effect of massive 
neutrinos on structure formation and on the CMB 
anisotropics. In section 4 we give an overview of 
recent cosmological neutrino mass limits and in 
Section 5 we discuss challenges for the future. 

2. THE EFFECT OF MASSIVE NEU- 
TRINOS ON STRUCTURE FORMA- 
TION 

The relic abundance of neutrinos in the Uni- 
verse today is straightforwardly found from the 
fact that they continue to follow the Fermi-Dirac 
distribution after freeze-out, and their tempera- 
ture is related to the CMB temperature Tcmb 
today by T, = (4/11)1/3Tcmb, giving 

_6C(3) 3 

where C(3) w 1.202, which gives n^, « 112 cm^^ 
at present. By now, massive neutrinos will have 
become non-relativistic, so that their present con- 
tribution to the mass density can be found by 
multiplying Ui, with the total mass of the neutri- 
nos m^^tot, giving 

" " 94eV' ^ > 

for Tcmb = 2.726 K. Several effects could modify 
this simple relation. If any of the neutrino chem- 
ical potentials were initially non-zero, or there 
were a sizable neutrino-antineutrino asymmetry, 
this would increase the energy density in neutri- 
nos and give an additional contribution to the 
relativistic energy density. However, from Big 
Bang Nucleosynthesis (BBN) one gets a very tight 
limit on the electron neutrino chemical potential, 
since the electron neutrino is directly involved 
in the processes that set the neutron-to-proton 
ratio. Also, within the standard three-neutrino 
framework one can extend this limit to the other 
flavours as well. Within the standard picture, 
equation should be accurate, and therefore 
any constraint on the cosmic mass density of neu- 
trinos should translate straightforwardly into a 
constraint on the total neutrino mass, according 



to equation If a fourth, light 'sterile' neu- 
trino exists, sterile-active oscillations would mod- 
ify this conclusion. Beacom et al. showed 
that extra couplings, not yet experimentally ex- 
cluded, of neutrinos may allow them to anni- 
hilate into light bosons at late times, and thus 
make a negligible contribution to the matter den- 
sity today. If so, equation ^ is not valid, and 
hence neutrino mass limits derived from large- 
scale structure do not apply. Wc shall assume 
that no such non-standard couplings of neutrinos 
exist. 

Finally, we assume that the neutrinos are 
nearly degenerate in mass. Current cosmologi- 
cal observations are sensitive to neutrino masses 
~ 1 eV or greater. Since the mass-square dif- 
ferences are small, the assumption of a degen- 
erate mass hierarchy is therefore justified. This 
is illustrated in Figure Q where we have plotted 
the mass eigenvalues mi, 7712,7713 as functions of 

ITly.tat ~ + ?772 + "^3 for AfTljl = 7 X 10^^ cV^ 

(solar) and Atti^j = 3 x 10"'^ eV^ (atmospheric), 
for the cases of a normal hierarchy (mi < m2 < 
TO3), and an inverted hierarchy (m3 < mi < m2). 
As seen in the Figure, for m^^tot > 0.4 eV the 
mass eigenvalues are essentially degenerate. 

Here we look at cosmological models with four 
components: baryons, cold dark matter, massive 
neutrinos, and a cosmological constant. Further- 
more, we restrict ourselves to adiabatic, linear 
perturbations. The basic physics is then fairly 
simple. Light, massive neutrinos can move un- 
hindered out of regions below a certain limiting 
length scale, and will therefore tend to damp a 
density perturbation at a rate which depends on 
their rms velocity. The presence of massive neu- 
trinos therefore introduces a new length scale, 
given by the size of the co-moving Jeans length 
when the neutrinos became non-relativistic. In 
terms of the comoving wavenumber, this is given 
by 

1 /2 

= 0.026 (^) ni(^hMpc~\ (3) 

for three equal-mass neutrinos, each with mass 
m^. The growth of Fourier modes with k > 
knr will be suppressed because of neutrino free- 
streaming. The free-streaming scale varies with 
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Figure 1. Neutrino mass eigenvalues as functions 
of TTiv^tat for the cases of normal (top panel) and 
inverted (bottom panel) hierarchies. The vertical 
line marked 'oscillations' is the lower limit derived 
from the measured mass-squared differences. The 
vertical line marked 'WMAP+SDSS' is the re- 
cent hmit derived in [SZ| from WMAP+SDSS, the 
vertical line marked '2dFGRS' is the limit from 
the 2dFGRS derived in "53 , and the line marked 
'■^H/3' is the upper limit from •^H /3 decay. 



the cosmological epoch, and the scale and time 
dependence of the power spectrum cannot be sep- 
arated, in contrast to the situation for models 
with cold dark matter only. 

The power spectrum of the matter fluctuations 
can be written as 

p„,(fc,z) -P,(fc)r2(fc,;2), (4) 

where T{k,z) is the 'transfer function', P^,{k) is 
the primordial spectrum of matter fluctuations, 
commonly assumed to be a simple power law 
P, (fc) = Ak", where A is the amplitude and the 
spectral index n is close to 1. It is also common to 
define power spectra for each component, see jl7| 
for a discussion. Note that the transfer functions 
and power spectra are independent of the value 
of the cosmological constant as long as it does 
not shift the epoch of matter-radiation equality 
significantly. 

The transfer function is found by solving the 
coupled fluid and Boltzmann equations for the 
various components. This can be done using one 
of the pubhcly available codes, e.g. CMBFAST 
HH] or CAMB [inj. In Figure|21we show the trans- 
fer functions for models with = 0.3, Qh = 
0.04, h — 0.7 held constant, but with varying 
neutrino mass m^. One can clearly see that the 
small-scale suppression of power becomes more 
pronounced as the neutrino fraction = fl^/il^ 
increases. 

The effect is also seen in the power spectrum, 
as shown in Figure |3| (top) . Note that the power 
spectra shown in the Figure have been convolved 
with the 2dFGRS window function, as described 
in |2()j . Furthermore, we have taken the possible 
bias of the distribution of galaxies with respect 
to that of the dark matter into account by leav- 
ing the overall amplitude of each power spectrum 
as a free parameter to be fitted to the 2dFGRS 
power spectrum data (the vertical bars in the Fig- 
ure). For a discussion of bias in the context of 
neutrino mass limits, see |21| . Because the errors 
on the data points are smaller at small scales, 
these points are given most weight in the fitting, 
and hence the power spectra in the Figure actu- 
ally deviate more and more from each other on 
large scales as increases. One can see from 
the Figure that a neutrino mass of riii, ~ 0.5 eV 
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Figure 2. Ratio of the transfer functions (at z = 
0) for various values of il^, to the one for fl^ — 0. 
The other parameters are fixed at ilni — 0.3, flh — 
0.04, h = 0.7. The solid line is for = 0.1 eV, 
the dashed line is for rrii, — 0.3 eV, the long- 
dashed line is for = 0.5 eV, and the dot- 
dashed line corresponds to nii, = 2 eV. 



or larger is in conflict with the data. The sup- 
pression of the power spectrum on small scales is 
roughly proportional to fi,: 



AlUk) 
Pm(fc) 



(5) 



This result can be derived from the equation of 
linear growth of density perturbations and the 
fact that only a fraction (1 — /^) of the matter 
can cluster when massive neutrinos are present 



3. CONSTRAINTS FROM THE CMB 
ALONE 

Neutrino masses also give rise to effects in the 
CMB power spectrum. If their masses are smaller 
than the temperature at recombination ^ 0.3 eV, 
their effect is very similar to that of massless neu- 




Figure 3. Top Figure: Power spectra for rrii, = 
(full hue), nil, = 0.1 (dotted hue), nii, = 0.3 
(dashed line), rrii, — 0.5 (long-dashed line), and 
nil, = 3 eV (dot-dashed line). The other param- 
eters are fixed at fim — 0.3, Ob = 0.04, h — 0.7. 
The vertical bars are the 2dFGRS power spec- 
trum data points. Bottom Figure: CMB power 
spectra for = (full line), rui, — 0.1 (dotted 
line), ni[y = 0.3 (dashed line), nii, — 0.5 (long- 
dashed line), and = 3 eV (dot-dashed line). 
The other parameters are fixed at Om = 0.3, 
Ob — 0.04, h = 0.7. The vertical bars are the 
WMAP power spectrum data points. 



5 



trinos IT . For slightly larger masses, there is an 
enhancement of the acoustic peaks with respect to 
the massless case, as shown in Figure |31 (bottom). 
While there is some sensitivity to the neutrino 
mass, note that all other parameters have been 
fixed in Figure O (bottom) . Analytic considera- 
tions by 24 provide insight into the effect of the 
neutrinos on the CMB. There are severe degen- 
eracies between rrii, and other parameters like n 
and nbh"^. The full analysis of the WMAP data 
alone in [26) gave no upper limit on m^. On the 
other hand |23] have claimed an upper limit of 2.2 
eV from CMB alone, in contrast with the conclu- 
sions of pri and [2S1- The differences might be 
due to the assumed priors and the marginalisation 
procedures over other cosmological parameters. 

Future CMB missions like Planck will provide 
high-resolution maps of the CMB temperature 
and polarization anisotropics. Gravitational lens- 
ing of these maps causes distortions, and Kapling- 
hat, Knox & Song |55] have shown that this ef- 
fect can be used to obtain very stringent lim- 
its on neutrino masses from the CMB alone. 
For Planck, they predict a sensitivity down to 
0.15 eV, whereas a future experiment with higher 
resolution and sensitivity can possibly reach the 
lower bound ~ 0.06 eV set by the neutrino oscil- 
lation experiments. 

4. RECENT COSMOLOGICAL NEU- 
TRINO MASS LIMITS 

The connection between neutrino masses and 
cosmic structure formation was realized early, but 
for a long time cosmologists were mostly inter- 
ested in neutrino masses in the ~ 10 eV range, 
since then they would be massive enough to make 
up all of the dark matter. The downfall of the 
top-down scenario of structure formation, and the 
fact that no evidence for neutrino masses existed 
before Super-Kamiokande detected oscillations of 
atmospheric neutrinos in 1998, makes it under- 
standable that there was very little continuous 
interest in this sub-field. However, the detection 
of neutrino oscillations showed that neutrinos in- 
deed have a mass. In an important paper Hu, 
Eisenstein & Tegmark showed that one could 
obtain useful upper limits on neutrino masses 



from a galaxy redshift survey of the size and qual- 
ity of the Sloan Digital Sky Survey (SDSS). 

Going down Table ^ one notes a marked im- 
provement in the constraints after the 2dFGRS 
power spectrum became available. After WMAP, 
there is a further tendency towards stronger up- 
per limits, reflecting the dual role of the CMB 
and large-scale structure in constraining neutrino 
masses: the matter power spectrum is most sensi- 
tive to neutrino masses, but one needs good con- 
straints on the other relevant cosmological pa- 
rameters to break degeneracies in order to ob- 
tain low upper limits. The limit will depend on 
the datasets and priors used in the analysis, but 
it seems like we are now converging to the pre- 
cision envisaged in The latest limit from 
\li7\ uses galaxy-galaxy lensing to extract infor- 
mation about the linear bias parameter in the 
SDSS, making a direct association between the 
galaxy and matter power spectra, and hence get 
a stronger constraint on the neutrino mass than 
would have been possible using just the shape of 
the galaxy power spectrum. 

Direct probes of the total matter distribution 
avoid the issue of bias and are therefore ideally 
suited for providing limits on the neutrino masses. 
Several ideas for how this can be done exist. In 
|29) the normalization of the matter power spec- 
trum on large scales derived from COBE was 
combined with constraints on erg (defined as the 
rms mass fluctuation in 8h~^Mpc radius sphere), 
they obtained a 95% confidence) from cluster 
abundances and a constraint m^^tot < 2.7 eV 
obtained, although with a fairly restricted pa- 
rameter space. However, as is probably one of 
the most debated numbers in cosmology at the 
moment |38| . and a better understanding of sys- 
tematic uncertainties connected with the various 
methods for extracting it from observations is 
needed before this method can provide useful con- 
straints. The potential of this method to push the 
value of the mass limit down also depends on the 
actual value of as- the higher as turns out to 
be, the less room there will be for massive neu- 
trinos. As an illustration we show in Figure 01 the 
value of as as a function of varying fl^, with the 
remaining cosmological parameters fixed at their 
'concordance' values. For a given value of m^^ 



6 



T? ATA T* ATI A A 




T SS 


V./ Liit;± 


limit 


lAAl 1 




Lya 


COBE norm., 
h = 0.72 ± 0.08, 

n cr£?oO 47 
(78 = 0.56ii„j ' 


5.5 eV 


ESI 






Jl,^ < 4 

" " 111 ^ W . , 

Qbh^ = 0.015, 
ft. < 0.8, n = 1.0 


2.7 eV 


EH 


pre-WMAP 


PSCz, Lya 




4.2 eV 


EH 


None 


2dFGRS 


BBN, SNIa, 
HST, n= 1.0 + 0.1 


2.2 eV 


EH 


pre-WMAP 


2dFGRS 




2.5 eV 


Ei] 


pre-WMAP 


2dFGRS 


SNIa, BBN 


0.9 eV 


m 


WMAP+CBI+ACBAR 


2dFGRS 


Lya 


0.71 cV 


Hi 


WMAP-I-Wang comp. 


2dFGRS 


HST, SNIa 


1.01 eV 


Em 


WMAP+CBI+ACBAR 


2dFGRS 


X-ray 


0.56t°io eV 


ESI 


WMAP 


SDSS 




1.7 eV 


ESI 


WMAP 


2dFGRS+SDSS 




0.75 eV 


ESI 


WMAP+ACBAR 


2dFGRS+SDSS 




1.0 eV 


EH 


WMAP 


SDSS 


bias 


0.54 eV 



one fits the corresponding CMB power spectrum 
to the data. This in turn leads to a best-fit ampli- 
tude and a prediction for erg for the given value of 
nil, ■ If one then has an independent measurement 
of (78, one can infer the value of m^. In Figure 
0]the amplitude of the power spectrum has been 
fixed by fitting to the WMAP data. The claimed 
detection of a non-zero neutrino mass in j30| can 
be seen to be due to the use of the cluster X- 
ray luminosity function to constrain erg, giving 
CTg = 0.69 ± 0.04 for = 0.3 [SSI- If a value 
of 178 at the higher end of the results reported in 
the hterature is used instead, e.g. as = 0.9 for 
flni = 0.3 from 40j , one gets a very tight upper 
limit on m^, but no detection of > 0. It is 
clearly important that systematic issues related 
to the various methods of obtaining as are set- 
tled. The evolution of cluster abundance with 
redshift may provide further constraints on neu- 
trino masses |41| . 

Direct probes of the mass distribution such as 
peculiar velocities and gravitational lensing are 
also potentially important for setting constraints 
on the neutrino mass. Deep and wide weak lens- 
ing surveys will in the future make it possible to 
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do weak lensing tomography of the matter den- 
sity field ^2 43 . By binning the galaxies in a 
deep and wide survey in redshift, one can probe 
the evolution of the gravitational potential. How- 
ever, because massive neutrinos and dark energy 
have similar effects on this evolution, complemen- 
tary information is required in order to break this 
degeneracy. Several studies of the potential of 
lensing tomography to constrain cosmological pa- 
rameters, in particular dark energy and neutrino 
masses, have been carried out, see e.g. for 
an overview. Even when taking the uncertain- 
ties in the properties of dark energy into account, 
the combination of weak lensing tomography and 
high-precision CMB experiments may be sensi- 
tive to neutrino masses below to lower bound of 
0.06 eV on the sum of the neutrino masses set by 
the current oscillation data j44| . 

5. DISCUSSION 

The dramatic increase in amount and quality of 
CMB and large-scale structure data we have seen 
in cosmology in the last few years have made it 
possible to derive fairly stringent limits on the 
neutrino mass scale. With the WMAP and SDSS 
data, the upper limit has been pushed down to 
~ 1 eV for the total mass, assuming three massive 
neutrino species. 

One point to bear in mind is that all these 
limits assume the 'concordance' ACDM model 
with adiabatic, scale-free primordial fluctuations. 
While the wealth of cosmological data strongly 
indicate that this is the correct basic picture, one 
should keep in mind that cosmological neutrino 
mass limits are mo del- dependent, and that there 
might still be surprises. As the suppression of 
the power spectrum depends on the ratio 

found that the out-of-fashion Mixed Dark 
Matter (MDM) model, with = 0.2, = 1 
and no cosmological constant, fits the 2dFGRS 
power spectrum well, but only for a Hubble con- 
stant Hq < 50 kms"""^ Mpc~^. A similar conclu- 
sion was reached in ^S] , and they also found that 
the CMB power spectrum could be fitted well by 
the same MDM model if one allows features in 
the primordial power spectrum. Another conse- 
quence of this is that excluding low values of the 



Hubble constant, e.g. with the HST Key Project, 
is important in order to get a strong upper limit 
on the neutrino masses. 

If the future observations live up to their 
promise, the prospects for pushing the cosmolog- 
ical neutrino mass limit down towards 0.1 eV are 
good. Then, as pointed out in "46" , one may even 
start to see effects of the different mass hierar- 
chies (normal or inverted), and thus one should 
take this into account when calculating CMB and 
matter power spectra. For example, with a non- 
degenerate mass hierarchy one will get more than 
one free-streaming scale, and this will leave an im- 
print on the matter power spectrum. The coming 
years will see further comparison between the ef- 
fective neutrino mass in Tritum beta decay, the 
effective Majorana neutrino mass in neutrinoless 
double beta decay and the sum of neutrino mass 
from Cosmology (J1Z|). It would be a great tri- 
umph for cosmology if the neutrino mass hierar- 
chy were finally revealed by the distribution of 
large-scale structures in the Universe. 
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